Review of self-potential methods in hydrogeophysics

The self-potential (SP) method is a passive geophysical method based on the natural occurrence of electrical fields on the Earth's surface. Combined with other geophysical methods, SP surveys are especially useful for localizing and quantifying groundwater flows and pollutant plume spreading, and estimating pertinent hydraulic properties of aquifers (water table, hydraulic conductivity). Laboratory experiments have shown that the involved coupling coefficients mainly depend on the fluid chemistry, conductivity and pH, and on the soil or rock properties. The interpretation of SP observations can be done qualitatively, for instance, by correlation of SP gradients with water fluxes (through electrokinetics) or salt fluxes (through electro-diffusion). In recent years, the interpretation has been improved with the help of modelling or/and inversion of the Poisson equation and endeavours to estimate hydraulic parameters by means of the intensity of electric current sources caused by underground flows.     

Assessing the Effect of Data Imports on the Completeness of OpenStreetMap – A United States Case Study

The assessment of OpenStreetMap (OSM) data quality has become an interdisciplinary research area over the recent years. The question of whether the OSM road network should be updated through periodic data imports from public domain data, or whether the currency of OSM data should rather rely on more traditional data collection efforts by active contributors, has led to perpetual debates within the OSM community. A US Census TIGER/Line 2005 import into OSM was accomplished in early 2008, which generated a road network foundation for the active community members in the US. In this study we perform a longitudinal analysis of road data for the US by comparing the development of OSM and TIGER/Line data since the initial TIGER/Line import. The analysis is performed for the 50 US states and the District of Columbia, and 70 Urbanized Areas. In almost all tested states and Urbanized Areas, OSM misses roads for motorized traffic when compared with TIGER/Line street data, while significant contributions could be observed in pedestrian related network data in OSM compared with corresponding TIGER/Line data. We conclude that the quality of OSM road data could be improved through new OSM editor tools allowing contributors to trace current TIGER/Line data.     

Tsunami hazard and exposure on the global scale

In the aftermath of the 2004 Indian Ocean tsunami, a large increase in the activity of tsunami hazard and risk mapping is observed. Most of these are site-specific studies with detailed modelling of the run-up locally. However, fewer studies exist on the regional and global scale. Therefore, tsunamis have been omitted in previous global studies comparing different natural hazards. Here, we present a first global tsunami hazard and population exposure study. A key topic is the development of a simple and robust method for obtaining reasonable estimates of the maximum water level during tsunami inundation. This method is mainly based on plane wave linear hydrostatic transect simulations, and validation against results from a standard run-up model is given. The global hazard study is scenario based, focusing on tsunamis caused by megathrust earthquakes only, as the largest events will often contribute more to the risk than the smaller events. Tsunamis caused by non-seismic sources are omitted. Hazard maps are implemented by conducting a number of tsunami scenario simulations supplemented with findings from literature. The maps are further used to quantify the number of people exposed to tsunamis using the Landscan population data set. Because of the large geographical extents, quantifying the tsunami hazard assessment is focusing on overall trends.     

Southern Hemisphere extra-tropical forcing: a new paradigm for El Ni?o-Southern Oscillation

The main goal of this paper is to shed additional light on the reciprocal dynamical linkages between mid-latitude Southern Hemisphere climate and the El Ni?o-Southern Oscillation (ENSO) signal. While our analysis confirms that ENSO is a dominant source of interannual variability in the Southern Hemisphere, it is also suggested here that subtropical dipole variability in both the Southern Indian and Atlantic Oceans triggered by Southern Hemisphere mid-latitude variability may also provide a controlling influence on ENSO in the equatorial Pacific. This subtropical forcing operates through various coupled air?Csea feedbacks involving the propagation of subtropical sea surface temperature (SST) anomalies into the deep tropics of the Atlantic and Indian Oceans from boreal winter to boreal spring and a subsequent dynamical atmospheric response to these SST anomalies linking the three tropical basins at the beginning of the boreal spring. This atmospheric response is characterized by a significant weakening of the equatorial Atlantic and Indian Inter-Tropical Convergence Zone (ITCZ). This weakened ITCZ forces an equatorial ??cold Kelvin wave?? response in the middle to upper troposphere that extends eastward from the heat sink regions into the western Pacific. By modulating the vertical temperature gradient and the stability of the atmosphere over the equatorial western Pacific Ocean, this Kelvin wave response promotes persistent zonal wind and convective anomalies over the western equatorial Pacific, which may trigger El Ni?o onset at the end of the boreal winter. These different processes explain why South Atlantic and Indian subtropical dipole time series indices are highly significant precursors of the Ni?o34 SST index several months in advance before the El Ni?o onset in the equatorial Pacific. This study illustrates that the atmospheric internal variability in the mid-latitudes of the Southern Hemisphere may significantly influence ENSO variability. However, this surprising relationship is observed only during recent decades, after the so-called 1976/1977 climate regime shift, suggesting a possible linkage with global warming or decadal fluctuations of the climate system.     

Extension tardi-orogénique et formation des bassins intracontinentaux: le bassin stéphanien des Cévennes

Résumé

Le bassin carbonifère des Cévennes, voisin du décrochement sénestre de Villefort, a été étudié en intégrant les données sédimentologiques, les données structurales du socle et du remplissage, ainsi que les données de la pétrologie métamorphique du socle. Le remplissage du bassin est contrôlé par des failles subméridiennes induisant une subsidence localisée en début d’ouverture. Cette subsidence devient plus régionale en fin d’histoire. Un grenat mangane- sifère est associé aux minéraux du faciès schistes verts, des blastes d’andalousite recoupent l’ensemble. Cette association minéralogique indique un gradient de température de 50°C/km. Ce gradient est semblable au gradient estimé dans le bassin sur la matière organique. La schistosité porte trois familles de linéation : 1- une famille à biotite-chlorite-quartz, orientée N45 ± 20° et associée à la mise en place des nappes cévenoles; 2- une famille à minéraux phylliteux dont l’orientation est située autour de N90 et associée à l’ouverture du bassin; 3- une famille à minéraux phylliteux orientée NO horizontale sur des plans de schistosité redressés près de la faille de Villefort. La schistosité régionale est affectée par des plis asymétriques, des bandes de cisaillement et des zones cata- clasées, le tout ayant une cinématique en faille normale vers l’Est. Les relations entre ces différentes structures suggèrent un continuum de déformation depuis des niveaux ductiles jusqu’à des niveaux fragiles. Le remplissage du bassin est affecté par des décrochements parallèles à la faille de Ville- fort. Le bassin des Cévennes s’est donc ouvert dans un contexte extensif est-ouest sur une croûte épaissie. La région a ensuite été reprise par des décrochements nord-sud dont celui de Villefort.     

Retention behaviour of humic substances onto mineral surfaces and consequences upon thorium (IV) mobility: case of iron oxides

The influence of humic acids (HA) on the retention of Th(IV) onto the surface of Fe oxides has been investigated. Thorium (IV) is here considered as an analogue of tetravalent radionuclides (Zr, Hf, U, Np, Pu). Hematite (α–Fe₂O₃) and gœthite (α–FeOOH) have been chosen as model phases because of their widespread presence in the earth's crust. Retention experiments were performed by batch procedure at constant ionic strength and varying pH in the range 2–10; the HA concentration varied in the range 1–100 mg l⁻¹. Compared with the former data obtained on the retention of Th(IV) with the same hematite sample free of HA, the study of a ternary system (i.e. hematite–HA–Th), clearly shows the influence of HA on the retention mechanisms of Th(IV) even in neutral to slightly alkaline medium. In particular, it is illustrated that the retention of Th(IV) onto Fe oxides closely depends upon the ratio between humic and surface (Fe oxide) sites. This ratio has to be greater than a critical value in order to impede the sorption of Th(IV).     

Receiver function method in reflection seismology

The receiver function method was originally developed to analyse earthquake data recorded by multicomponent (3C) sensors and consists in deconvolving the horizontal component by the vertical component. The deconvolution process removes travel path effects from the source to the base of the target as well as the earthquake source signature. In addition, it provides the possibility of separating the emergent P and PS waves based on adaptive subtraction between recorded components if plane waves of constant ray parameters are considered. The resulting receiver function signal is the local PS wave's impulse response generated at impedance contrasts below the 3C receiver.We propose to adapt this technique to the wide‐angle multi‐component reflection acquisition geometry. We focus on the simplest case of land data reflection acquisition. Our adapted version of the receiver function approach consists in a multi‐step procedure that first removes the P wavefield recorded on the horizontal component and next removes the source signature. The separation step is performed in the τ?p domain while the source designature can be achieved in either the τ?p or the t?x domain. Our technique does not require any a priori knowledge of the subsurface. The resulting receiver function is a pure PS‐wave reflectivity response, which can be used for amplitude versus slowness or offset analysis. Stack of the receiver function leads to a high‐quality S wave image.     

Contrasting regional discharge evolutions in the Amazon basin (1974–2004)

Former hydrological studies in the Amazon Basin generally describe annual discharge variability on the main stem. However, the downstream Amazon River only represents the mean state of the Amazonian hydrological system. This study therefore uses a new data set including daily discharge in 18 sub-basins to analyze the variability of regional extremes in the Amazon basin, after recalling the diversity of the hydrological annual cycles within the Amazon basin. Several statistical tests are applied in order to detect trends and breaks in the time series. We show that during the 1974–2004 period, the stability of the mean discharge on the main stem in Óbidos is explained by opposite regional features that principally involve Andean rivers: a decrease in the low stage runoff, particularly important in the southern regions, and an increase in the high stage runoff in the northwestern region. Both features are observed from the beginning of the nineties. These features are also observed in smaller meridian sub-basins in Peru and Bolivia. Moreover we show that the changes in discharge extremes are related to the regional pluriannual rainfall variability and the associated atmospheric circulation as well as to tropical large-scale climatic indicators.     

An external–internal mode coupling for a 3D hydrodynamical model for applications at regional scale (MARS)

This paper presents a 3D model in sigma coordinates. Although the principles it is based on have been established for some time, some original aspects for this type of 3D mode splitting model are presented here. The model was designed to simulate flows in various coastal areas from the regional scale down to the inshore scale of small bays or estuaries where circulation is generally driven by a mix of processes. The processes to be modeled enable simplifications of the Navier–Stokes equations on the classic Boussinesq and hydrostatic hypotheses. These equations are transformed within a sigma framework to make free surface processing easier. The main point of our demonstration focuses on the original aspect of the coupling between barotropic and baroclinic modes especially designed for ADI. It explains how full consistency of the transport calculated within the 2D and 3D equation sets was obtained. Lastly, we describe the physical processes simulated on a realistic configuration at a regional scale in the Bay of Biscay.